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Recent obserbations by Fermi Gamma-Ray Space Telescope of 7-ray pulsars have revealed further details of 
the structure of the emission region. We investigate the emission region for the multi-wavelength light curve 
using outer gap model. We assume that 7-ray and non-thermal X-ray photons are emitted from a particle 
acceleration region in the outer magnetosphere, and UV/optical photons originate above that region. We also 
assume that 7-rays are radiated only by outwardly moving particles, whereas the other photons are produced 
by particles moving inward and outward. We parametrize the altitude of the emission region. We find that the 
outer gap model can explain the multi-wavelength pulse behavior. From observational fitting, we also find a 
general tendency for the altitude of the 7-ray emission region to depend on the inclination angle. In particular, 
the emission region for low inclination angle is required to be located in very low altitude, which corresponds to 
the inner region within the last-open field line of rotating dipole in vacuum. This model suggests a modification 
of statistics about observed 7-ray pulsars. 



1. Introduction 

Pulsars emit over a wide range of energies from 
radio to 7-ray. Recent observations by the Fermi 
Gamma-Ray Space Telescope of more than eighty pul- 
sars i^] have revealed further details of the structure 
of the emission region. The light curve in the 7-ray 
band is an important tool for probing the particle ac- 
celeration and dissipation processes in the pulsar mag- 
netosphere. The 7-ray emission region has therefore 
been explored by comparing theoretical models with 
the observed light curve(e.g., [HI [13, ) • The pulsed 
emission is also detected in other energy bands (X-ray, 
ultraviolet, optical and radio) for some sources. The 
spectral features are non-thermal except for the soft 
X-ray range, and the light curves from a single ob- 
ject are, in general, different from one energy band to 
another. A complete understanding of light curve be- 
havior in multi- wavelength bands can provide valuable 
information about the particle acceleration region. 

Possible origins of non-thermal pulsed emissions 
have been considered in the polar cap ^ , slot gap [l3| , 
and outer gap [^ models. Recent Fermi observations 
[l[ rules out the near-surface emission proposed in po- 
lar cap cascade models, which would exhibit a much 
sharp spectral cutoff due to magnetic pair-production 
attenuation. Thus, pulsed 7-ray emission originates in 
the outer magnetosphere, as considered in the outer 
gap model. 

Here, we investigate the emission regions of several 
pulsars by fitting the simplified model of Takata et al. 
(2008, hereafter TCS08) [H to the observed muhi- 
wavelength light curves. In this model, we have to 
specify the locations of the upper and lower bound- 
aries of the gap region where the non-corotation po- 
tential is zero. Therefore, we explicitly introduce the 
altitude of the gap region as a parameter, in order 
to fit the observational data easily. The light curves 
also depend on the dipole inclination and viewing an- 
gles. In our method, such parameters are eliminated 
by other observational data, and only the altitude is 



changed for the fitting. In the most studies, the lower 
boundary of the emission region is chosen as the sur- 
face of the last-open field lines of the rotating dipole 
(e.g., [14]). In this study, however, the altitude is 
allowed to be in a wide range in order to explore the 
possible deviation of magnetic field-line structure from 
that of a rotating dipole in vacuum. 



2. Model 

The numerical method for fitting the light curve is 
well described by Romani & Watters (2010) [lT| and 
Bai & Spitkovsky (2010a) [3]. We follow this method 
except one modification that we explicitly introduce 
the altitude of the emission region as an additional 
parameter. We assume that magnetic field structure 
is approximately described by a rotating dipole and 
that radiation direction aligns with magnetic field in 
a frame rotating with angular velocity in which the 
electric field vanishes. 

A certain mechanism is needed to fix the lower 
boundary of the particle acceleration region. In most 
works, including TCS08, the lower boundary is chosen 
as the surface of the last-open field lines of a rotating 
dipole in a vacuum. In the outer gap model, if particle 
acceleration occurs in an open zone, the curvature ra- 
diation from the accelerated particles forms a narrow 
cone along the magnetic field lines in a frame rotating 
with angular velocity. These 7-ray photons are con- 
verted by colliding X-ray photons to pairs, which 
tend to screen the accelerating electric field. How- 
ever, there is no supply of pairs on the last-open field 
lines and hence no screening of the electric field, since 
the 7-ray photons are emitted only toward higher al- 
titudes above the last-open field lines [5]. The 'real' 
last-open field lines may be different from ones in a 
vacuum. We therefore take into account this possi- 
ble deviation of the boundary. We assume that dipole 
magnetic field is an approximation within the light 
cylinder and use rotating diople field as the global 
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magnetic field structure. Even if the overall structure 
is not different so much, critical value between open 
and closed field lines is very sensitive to the boundary 
value at the surface. Thus we introduce a parameter, 
altitude of the emission region as a correction factor in 
order to take into account the deviation of boundary 
from the vacuum field. In our model this parameter 
specifies the range of the emission region which is lo- 
cated above or below the last-open field lines within 
the light cylinder radius. Each different field line orig- 
inating from the magnetic polar region is parameter- 
ized by magnetic colatitudes and azimuthal angles 
0in- Following Cheng et al. (2000) [6^, we define open 
volume coordinates on the polar cap, (rov, 0m) i where 
rov = 0^/e^-^{(j)^). The function 9^^° is the mag- 
netic colatitude of the conventional polar cap angle 
and generally depends on the magnetic azimuth 
The parameter rov corresponds to the altitude of the 
emission region: The last-open field lines of a rotat- 
ing dipole in a vacuum correspond to Tqv = 1, whereas 
those for higher altitudes have rov < 1 • The maximum 
value is chosen as rov = l.SB'^/^, which corresponds to 
the polar cap angle 0"^ ^ 1.36^/^0Pf'°, obtained in the 
force-free limit by Contopoulos et al. (1999) 0- 

We assume that the radiation of different energy 
bands is emitted from different field lines character- 
ized by altitude. The field line relevant to the 7-ray 
and X-ray is approximated as being the same one. The 
direction of the emission is tangential to the lines, and 
inward and outward directions are possible. Both lo- 
cation and direction affect the light curve profile of 
the energy bands. Following the model by TCS08, 
the 7-ray radiation above 100 MeV is emitted by par- 
ticles moving in an outward direction, whereas radia- 
tion at lower energy bands is emitted by those mov- 
ing in both outward and inward directions. We use 
two conditions to constrain the emission region. First 
condition is the radial extension of the emission re- 
gion. The outward emission is restricted to radial dis- 
tances rn < i?LCj and the inward one is restricted to 
Ts < r < min(3rn, i?Lc)- The outer boundary 3rn for 
inward emission comes from the results of dynamic 
model (TCS08), in which very few ingoing pairs are 
produced beyond the radial distance r > 3rn. Second 
condition is the azimuthal extension of the emission 
region. We use the magnetic azimuthal angle of the 
footprint of field line (i.e., the point where magnetic 
field line penetrates the neutron star surface) to char- 
acterize the field line for given rov. Radial distance 
to the null charge surface on the field lines signifi- 
cantly depends on the magnetic azimuthal angle. In 
the outer gap model, most of the pairs are created 
around the null surface (TCS08). We expect that the 
gap activity is related to the distance to the null sur- 
face. Although the current density should be deter- 
mined by global conditions, there is no study of the 
three-dimensional magnetosphere of an inclined rota- 
tor. Here, we assume that the field lines of both out- 



ward and inward emission are active only if the radial 
distance to null surface r,; is shorter than i?LC- Spa- 
tial distribution of the emissivity is approximated by 
the step function- type, but the peak positions weakly 
depend on the detailed emissivity distribution. 

We assume that the overall structure of the light 
curve comes not from the emissivity distribution, but 
from a bunch of many field lines in the observation, 
that is, caustics. Here, we focus on the peak phases of 
the light curve, so we adopt a simple, uniform emissiv- 
ity along all magnetic field lines, which is independent 
of both the magnetic azimuthal angle 0m and the al- 
titude rov. 

In the observed light curve, the reference phase 
= is assumed to be located at the radio emission 
peak maximum. However in the model light curve, 
the conventional reference phase = occurs when 
the magnetic axis, spin axis and Earth line of sight 
lie all in the same plane. These two reference phases 
do not agree with each other since it is generally as- 
sumed that radio emissions arise at non-zero altitude 
in most empirical studies. Following Romani & Wat- 
ters (2010) [ll|, we allow a shift by -0.1 < Scf) < 0.1 
in the model reference rotation phase. This degree of 
freedom does not significantly affect the determina- 
tion of the altitude parameter rov. For more details 
on this model see Kisaka & Kojima (2011) 

3. Results 

Here, we compare our model with pulse profiles ob- 
served at multiple wavelengths for seven pulsars. The 
sources are chosen using two criteria. One is that non- 
thermal pulses are detected in addition to the 7-ray 
and radio bands. Our concern is to explore whether 
or not the emission region for different energy bands is 
explained by the TCS08 model. The second criterion 
is that the geometrical parameters, inclination angle 
a and viewing angle ^ are observationally constrained 
by the relativistic Doppler-boosted X-ray pulsar wind 
nebula or radio polarization data. The torus fitting 
method constrains the viewing angle ^ only. A small 
allowed range of [a — ^| < 10° is assumed for samples 
in which only ^ is constrained due to the fact that ra- 
dio emission from the pulsar polar region is detected. 
The geometrical parameters for our selected pulsars 
are listed in Table 1. We use these values, although 
there are some uncertainties in them. The results are 
summarized in Fig. 1, which shows the intensity map 
for outward (upper panel) and inward (lower panel) 
emission as a function of the altitude of the emission 
region rov and rotational phase. This shows that the 
peak phases of seven pulsars emitting 7-ray and X-ray 
can be successfully fitted using the TCS08 outer gap 
model, in which both 7-rays and X-rays originate from 
the same magnetic field line characterized by an alti- 
tude rov ■ The parameter rov > 1 is needed in the light 
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TABLE 1 Geometrtical parameters 



Name 


a 






(degrees) 


(degrees) 


(1) 


(2) 


(3) 


J0835-4510 


72 


64 


J0659+1414 


29 


38 


J0205+6449 


78 


88 


J2229+6114 


55 


46 


J1420-6048 


30 


35 


J2021+3651 


75 


85 


J1057-5226 


75 


69 




Table I NOTES.-Col.(l):Pulsar name.Col.(2):The 
inclination angle. Col.(3):The viewing angle. 



Xoptoptx Topt. xopl X X 
(A) X I II (B) 4 4 (C) I I 




0.00.20.406081 00020406081 00204060 81 00020406081 



Figure 1: The intensity maps for seven pulsars. Upper 
left to right: PSRs J0835-4510 (A), J0659+1414 (B) and 
J0205+6449 (C). Lower left to right: PSRs J2229+6114 
(D), J1420-6048 (E), J2021+3651 (F) and J1057-5226 
(G). In each sample, upper panel is outward emission 
and lower is inward emission. The blue horizontal arrows 
show best fit values of rov for 7-ray and X-ray emission 
regions. The red vertical arrows show the phase of peaks. 
The red horizontal arrows in (E) and (G) show the phase 
range of broad peaks. See Kisaka & Kojima (2011) .9]. 



curve fitting for some sources. Moreover, the inclusion 
of inward emission for X-rays causes a variety of pulse 
profiles in both bands. The parameter Tov could not 
be determined solely using 7-ray data for a single 7- 
ray peak pulsar. However, by considering the X-ray 
light curve, the parameter is uniquely determined for 
PSRs J0659-M414, J2229-h6114 and J1420-6048. For 
details of each pulsars, see Kisaka & Kojima (2011) 
i. 

It is worthwhile to explore the general dependence 
of the altitude rov on other characteristics if any, al- 
though there may not be enough data for a proper 
statistical analysis. In Fig. 2, Tov is plotted as a func- 
tion of inclination angle a, spin-down luminosity Lgd, 
characteristic age Tc and surface dipole magnetic field 
Bs. We found that there is a significant correlation 
between and the inclination angle a only; the rela- 



Figure 2: The relation between rov and inclination angle 
(upper), spin-down luminosity (lower left), characteristic 
age (lower middle) and surface magnetic field (lower 
right). The altitudes corresponding to the separatrix 
layer model are shown as purple downward and blue 
upward triangles in the upper panel. The two lines are 
linear fitting lines for the separatrix layer model. The 
light green curve in the lower left panel shows the 
relation (1 - rov) = (10^^ erg /L^dY^^. See Kisaka & 
Kojima (2011) 3. 



tions of rov with the other parameters are very weak. 
This correlation suggests that the deviation from a 
vacuum rotating dipole field is large for small incli- 
nation angle. It is very interesting to compare this 
result with that in a force-free magnetosphere. Bai & 
Spitkovsky (2010b) [4] proposed that the separatrix 
layer at an altitude of 0.90-0.95 times the height of 
the last-open field line is relevant to emissions in a 
three-dimensional inclined force-free magnetosphere. 
This altitude, which is not exactly symmetric with re- 
spect to the magnetic azimuthal angle but can 
be approximated by the value at (/)ni = 0, is plotted in 
Fig. 2 as purple downward and blue upward triangles. 
Two linear fitting lines are also shown. The altitude 
rov decreases with the inclination angle a in both our 
model and the separatrix layer model of a force-free 
magnetosphere. 

The thickness of the gap region, w, is not known, 
but it is sometimes assumed to decrease with the spin- 
down luminosity Lgd [HI, [3 • We have w = 1 — rov , if 
the lower boundary of the gap is fixed as the last-open 
field line in the vacuum dipole field. This assumption 
is tested in the lower left panel of Fig. 2, in which the 
relation (1 - rov) ~ (^sd/lO^^erg s~^)"^/^ is plotted 
as a light green curve. (The curve is not fitted to 
the data points.) This suggests that the assumption 
of maximum altitude, rov = 1-0, is not a good one. 
This discovery affects expected number of the 7-ray 
pulsars in the observation. From geometrical reason, 
the pulsed emission by caustics is limited to a certain 
range between inclination and viewing angles. 

Romani & Walters (2010) [Hj showed the range of 
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30 60 90 



dipole field. The real last-open field lines lie inside 
those of vacuum dipole field, Tov < 1-0. This shift 
suggests that the lower boundary is very similar to 
that of a force-free magnetosphere. We find that the 
altitude of the emission region is correlated with incli- 
nation angle. This relationship is also very similar to 
that in a force-free magnetosphere. The lower bound- 
ary of emission region has been assumed to = 1 so 
far, but our model fits do not support it. This mod- 
ification of the boundary of the magnetosphere sug- 
gests that the pulsars with low inclination and view- 
ing angles are likely to be detectable. Thus the ex- 
pected number in the future observation in the previ- 
ous works is underestimated for the sources with low 
inclination and viewing angles. 



Figure 3: The observable range for 7-ray pulsars in the 
a-f plane for the outer gap model. Black solid curve 
shows the boundary of observable pulsars using linear 
fitting line for the separatrix layer model. Red dashed, 
blue dash-dotted and purple dotted curves show the 
boundary with rov =0.95, 0.90 and 0.70, respectively. 
Light-green curcles show the pulsars in Table 1. See 
Kisaka & Kojima (2011) [1]. 



observable pulsars with =0.95, 0.90 and 0.70 for 
outer gap model in their Fig. 16. We recalculate it 
and show the result in Fig. 3. The observable range 
of viewing angle ^ is below the curves. Our finding 
in Fig. 2 is that is a function of the inclination 
angle, which is similar to that of the separatrix layer 
model. We also show the observable range by the em- 
pirical relation obtained in Fig. 2 as black solid line, 
for which the altitude is chosen as 0.925 times the 
height of the last-open field line in force-free magne- 
tosphere. The figure shows that sources with low in- 
clination and viewing angles become observable. For 
example, pulsar with the inclination angle a = 30° 
can be detected for ^ > 60° for Tov =0.95, but for 
^ > 30°. Thus expected number increases approxi- 
mately twice for sources with the low inclination and 
viewing angles. 



4. SUMMARY 

We have calculated the light curves of emissions us- 
ing the TCS08 outer gap model and compared them 
with observed multi- wavelength light curves. We find 
that the model can successfully explain the peak po- 
sitions of multi-wavelength light curves. In order to 
determine the altitude of the emission region, the ob- 
served X-ray light curve is important, especially when 
there is a single peak in the 7-ray light curve. 

The best-fit values of the altitude of the emission 
region for PSRs J0659-M414 and J1420-6048, suggest 
a deviation from the last-open field lines of a vacuum 
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